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An analysis is carried out within mSUGRA of the estimated number of events originating from 
upward moving ultra-high energy neutralinos passing through Earth's crust that could be detected 
by the Extreme Universe Space Observatory (EUSO). The analysis exploits a recently proposed 
technique that differentiates ultra-high energy neutralinos from ultra-high energy neutrinos using 
their different absorption lengths in the Earth's crust. It is shown that for the part of the parameter 
space, where the neutralino is mostly a Bino and with squark mass ~ 1 TeV, EUSO could see ultra- 
high energy neutralino events within mSUGRA models with essentially no background. In the energy 
range 10 9 GeV < E% < 10 11 GeV the unprecedented aperture of EUSO makes the telescope sensitive, 
after 3 yr of observation, to neutralino fluxes as low as d^/dE^ > LI x 10~ 6 (E-^/GeV) -1 ' 3 GeV -1 
cm -2 yr -1 sr _1 , at the 95% CL. Such a hard spectrum is characteristic of supermassive particles' 
iV-body hadronic decay. The case in which the flux of ultra-high energy neutralinos is produced via 
decay of metastable heavy (mx — 2 x I0 12 GeV) particles with uniform distribution throughout the 
universe, and primary decay mode into 5 quarks + 5 squarks, is analyzed in detail. The normalization 
of the ratio of the relics' density to their lifetime has been fixed so that the baryon flux produced 
in the supermassive particle decays contributes to about 1/3 of the events reported by the AG AS A 
Collaboration below 10 11 GeV, and hence the associated GeV 7-ray flux is in complete agreement 
with EGRET data. For this particular case, EUSO will collect between 4 and 5 neutralino events 
(with 0.3 of background) in m 3 yr of running. NASA's planned mission, the Orbiting Wide-angle 
Light-collectors (OWL), is also briefly discussed in this context. 

PACS numbers: 11.30.Pb, 96.40.Pq - NUB-3247/Th-04 



I. INTRODUCTION 



mSUGRA [l] and its extensions (generically called 
SUGRA models) are currently the leading candidates 
for physics beyond the standard model. These models 
contain a consistent mechanism for the breaking of su- 
persymmetry softly by gravity mediation. An attractive 
feature of these models is that with R parity conserva- 
tion the lightest neutralino is a possible candidate for 
cold dark matter in a significant part of the mSUGRA 
parameter space Further, over most of the parame- 
ter space the phenomenon of scaling occurs |3| so that 
the light neutralino is mostly the supersymmetric part- 
ner of the U(1)y gauge boson B^, i.e., it is mostly a 
U(1)y gaugino or a Bino 0,0. The parameter space of 
mSUGRA is characterized by the universal scalar mass, 
777.0 , the universal gaugino mass, mi , the universal trilin- 
ear coupling, A (all taken at the grand unification scale 
M G ~ 2 x 10 16 GeV), tan/3 =< H 2 > / < H x > where 
H 2 gives mass to the up quark, and Hi gives mass to the 
down quark and the lepton. In addition the model con- 
tains the Higgs mixing parameter /i which enters in the 
superpotential in the form \xH\Hi. The magnitude of /x 
is determined by the constraint of radiative electro-weak 
symmetry breaking in the theory while signal is arbitrary 
and must be constrained by experiment. 
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mSUGRA has been put to stringent test by the re- 
cent precision data from the satellite experiment, the 
Wilkinson Microwave Anisotropy Probe (WMAP) which 
imposes a narrow range for cold dark matter (CDM) 
so that H @ n C DMh 2 = 0-1126±g;gg|. The candi- 
dacy of neutralinos as the dark matter of the universe is 
based on relic densities surviving annihilation processes 
of non-relativistic particle. Detailed analyses show that 
mSUGRA allows for a small amount of the parameter 
space in agreement with WMAP observations Q. As 
a consequence, the applicability of mSUGRA demands 
that the contribution of other sources of CDM to the 
dark matter mix are negligible. In this paper, we will 
be interested in a flux of ultrarelativistic neutralinos re- 
sulting from decays of apopulation of CDM metastable 
superheavy particles 0, ||J . In concert with the previous 
statement, these particles should contribute negligibly to 
the dark matter density. 

The weak couplings of neutralinos imply an interaction 
length in air which is greater than the atmospheric depth, 
even at horizontal incidence. The interaction probability 
is then roughly uniform throughout the atmosphere. As 
with neutrinos, showers initiated by neutralino primaries 
can be distinguished from hadronic events by restricting 
the zenith angle space to near horizontal - this maximizes 
the probability to detect showers of weakly interacting 
primaries, while screening out the electromagnetic com- 
ponent of hadronic showers which are initiated high in 
the atmosphere. However, deeply developing neutralino 
cascades cannot be isolated from neutrino induced air 
showers. 

In this paper we show that the part of the parame- 
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ter space where the neutralino is mostly a Bino and the 
mass Trig of the first and second generation squarks is 
~ f TeV can lead to ultra-high energy neutralino sig- 
nals that may be see n by the Extreme Universe Space 
Observatory (EUSO) }wL UA\ . These two conditions are 
fully compatible with the WMAP constraint, and the 
neutralino as lightest supersymmetric particle. Further, 
with appropriate cuts the background events arising from 
ultra-high energy neutrinos are essentially negligible. We 
discuss now the details of the analysis. 

The problem of discriminating between neutralino and 
neutrino induced showers with space-based experiments 
has been examined recently ^3 • The method makes use 
of the Earth as a filter. Specifically, in the region of the 
mSUGRA parameter space under consideration the cross 
section for neutralino-nucleon interaction is smaller than 
that for neutrino-nucleon scattering processes. Thus, 
by restricting the angular bin for arrival of upward go- 
ing showers to a region where neutrinos are largely ab- 
sorbed during traversal in the Earth, it may be possible 
to obtain a clean signal In Ref. the discus- 

sion was presented in terms of neutralino-nucleon cross 
sections parameterized as a series of constant fractions 
of the neutrino-nucleon cross section. In this paper, we 
first calculate the neutralino cross section in the squark- 
resonance approximation. We then proceed to estimate 
the sensitivity of EUSO to ncutralino-induced air show- 
ers. The sensitivity will be characterized by a lower 
bound on the neutralino flux, which is then related to 
some particular models of X-particle decay. 



II. EUSO OUTLINE 

In view of a mission starting in 2006-2007, the Extreme 
Universe Space Observatory has been designed to 

observe the complex relativistic cascades induced by the 
incoming extra terrestrial radiation using sensors in the 
ultraviolet band (300 - 400 nm), a technique pioneered by 
the Fly's Eye experiment ^4| . The telescope will operate 
for more than 3 yr on board of the International Space 
Station. The eye will be equipped with wide-angle Fres- 
nel optics lens that provide a field-of-view of ±30°, at an 
orbit altitude of about 400 km. The monocular stand- 
alone configuration of the instrument will cover an area of 
w 1.6 x 10 5 km 2 , imaging an air target mass that exceeds 
10 12 tons. This corresponds to a water equivalent (w.e.) 
effective volume of ~ 2400 km 3 . Observations can only 
be done on clear moonless nights. Limitations associated 
with the cloud system and ultraviolet background sources 
result in an average 10% - 15% duty cycle 0. Hence, 
the incredibly large geometric aperture, A w 7500 km 3 
w.e. sr, is somewhat reduced. 

The fluorescence eye consists of several large light col- 
lectors (or telescopes) which image regions of the sky onto 
clusters of light sensing and amplification devices. The 
basic elements of a telescope are the diaphragm, which 
defines the telescope aperture, the spherical mirror that 



must be dimensioned to collect all the light entering the 
diaphragm in the acceptance angular range, and the cam- 
era which consists of an array of photomultiplier tubes 
(PMTs) positioned approximately on the mirror focal 
surface. The PMTs effectively pixelize the region of the 
sky covered by the telescope. The shower development is 
detected as a long, rather narrow sequence of hit PMTs. 
As the point-like image of the shower proceeds through 
an individual PMT, the signal rises, levels off, and falls 
again. 

The sensitivity of the detector depends primarily on 
the signal (S) to noise (N) ratio [l4j . 
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where (B) w 400 m~ 2 sr -1 ns _1 is the average pho- 
ton night glow background in the ultraviolet band [lRj . 
e « 20% is the quantum efficiency of the detector, Ki w 1 
is the transmission coefficient in the ozone layer, ki = 



-r/\ R 



is the transmission coefficient in the atmosphere, 



Ar ~ 23 km (at 400 nm) is the Rayleigh scattering length 
at sea level, R p is the distance of closest approach be- 
tween the shower and the detector, r w 15 km is the 
effective slant depth of the lower atmosphere, D — 2.5 m 
is the diameter of the mirror aperture, and d = 5 mm is 
the diameter of the PMT. The fluorescence trail is emit- 
ted isotropically with an intensity that is proportional to 
the number of charged particles in the shower, N e . The 
fluorescent signal is roughly N 7 ~ 4 photons/electron/m 
and N e m 0.625 (E/GeV). This translates into a 4cr en- 
ergy threshold of 

E = 2.7x 10 6 i? 3/2 exp(r/A_R) GeV . (2) 

In addition to the fluorescence process, the electrons pro- 
duce a large photon signal from Cerenkov radiation that 
is primarily beamed in the forward direction. The pro- 
duction of fluorescence light is less than 10 -4 of that in 
the Cerenkov cone [l6|. For Earth-penetrating showers 
emerging upward in the direction of the orbiting tele- 
scope, this Cerenkov signal extends the 4cr threshold of 
EUSO to the PeV energy band. 



III. CROSS SECTIONS 

The cross section for the resonant scattering x+q(q) 
q(q) all [HI 



a xN = V 2tt / q(x) — ^ — 

x i^ianr 2 ^^) 



(3) 



where p q and p^- are the 4-momenta of the quark and 
neutralino in the interaction, 
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where we have included the contributions from the left- 
handed and right handed squarks assuming they are de- 
generate. In the above the first sum runs over all quark 
flavors with parton distribution functions (pdf's) indi- 
cated by q(x). For small L-R mixing, and ignoring quark 
masses and small Higgsino mixings one has ftfty 



^= = (X 2 cos9 w - Xi sin 0^0 (J^ - sin 2 9 W 



where, using the lab relation pn-Px — MpjE x , we now 
have 
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We use the CTEQ6D pdf's H3, which in the energy 
region of interest for up, down, strange and charmed 
(quarks and anti quarks) can be parametrized as 
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and similar relations hold for the charm and the strange 
squarks. In Eqs. ||SJ, @), Q, and © ff2 is the weak 
SU(2) gauge coupling constant, 9w is the weak mixing 
angle, e = g 2 sin9w is the U(l) cm charge, and X\,X 2 
are the projection of x on t ne Bino and Wino super 
partners, respectively, and we have ignored small Hig- 
gsino components X 3 and X4 where Ylt=i \Xi\ 2 = 1. 
In the analysis we use renormalization group evolution 
of the soft parameters from the grand unification scale 
Mc to low energy to compute the neutralino mass ma- 
trix. This matrix depends on the sign of /i. The recent 
experiment on the anomalous magnetic moment of the 
muonjljj indicates a positive sign for fj, for the supersym- 
metric contribution|20| and thus in this analysis we have 
chosen a positive fj, sign. The diagonalization of the neu- 
tralino mass matrix determines the projections Xi,X 2 
which are computed in terms of the input mSUGRA pa- 
rameters. In the parameter space of interest, i.e., where 
rriq ~ 1 TeV, and the WMAP constraint is obeyed, the 
neutralino is essentially a Bino, and one has X\ ~ 1 while 
X 2 , X3, X4 are relatively small. 

Through use of the parton model relation 2p q .p x = 
2xPn-Px {Pn is the nucleoli momentum) and the 5- 
function, one finds the compact expression 
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Then, for 

VlTeV/ 1 



GeV VlTevJ ' [ ' 

the neutralino-nucleon interaction cross section becomes 
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We comment briefly on the relation of our calculation 
to a previous estimate j 2 21 . In contrast to that work, we 
do not characterize the SUSY parameter space accord- 
ing to the decay branching ratio T(q — » qg)/T(q — > q\). 
The total cross section, as given above, is independent of 
any branching to gluinos. Comparison of Eq. I|13f) with 
the cross section for the competing process X g ^ qq es- 
timated in [22| shows that the resonant cross section is 
about an order of magnitude greater. 



IV. EVENT RATES 

The most popular mechanism to date to produce ultra- 
high energy neutralinos is annihilation or decay of super 
massive (10 12 GeV < m x < 10 16 GeV) V-particles. To 
maintain an appreciable decay rate today, one has to tune 
the X lifetime to be longer (but not too much longer) 
than the age of the universe |8|, |9| , or else "store" short- 
lived V-particles in topological vestiges of early universe 
phase transitions |2^|. The cascade decay to cosmic ray 
particles is driven by the ratio of the volume density of 
the V-particle (nx — p c &x /mx) to its decay time (r x ). 
This is very model dependent, as neither the cosmic aver- 
age mass density contributed by the relics (Qx), nor r x is 
known with any degree of confidence (p c w 1.05 x 10~ 4 h 2 
GeV cm -3 , with h = Hubble constant in units of 100 km 
sec -1 Mpc" 1 ). Moreover, the internal mechanisms of the 
decay and the detailed dynamics of the first secondaries 
do depend on the exact nature of the particles. Conse- 
quently, no firm prediction on the expected flux of neu- 
tralinos can be made. However, if there are no new mass 
scales between MgusY ~ 1 TeV and mx , the squark and 
sleptons would behave like their corresponding supersym- 
metric partners, enabling one to infer from the "known" 
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evolution of quarks and leptons the gross features of the 
A-particle cascade: the quarks hadronize producing jets 
of hadrons containing mainly pions together with a 3% 
admixture of nucleons p4j |. The final spectrum is then 
dominated by gamma rays and neutrinos produced via 
pion decay. 

In light of the mounting evidence that ultra-high en- 
ergy cosmic rays are not gamma rays |25| . proton dom- 
inance at ultra-high energies is achieved by efficient ab- 
sorption of the dominant ultra-high energy photon flux 
on the universal and/or Galactic radio background. This 
results in a recycling of the photon energy down to the 
MeV - GeV region. Thus, since the baryon and photon 
components of the A-particle decay are correlated by the 
dynamics, the measurement of the GeV diffuse gamma 
ray flux can significantly constrain the cosmic ray pro- 
duction by A"-particles, integrated over redshift [26j. In 
this direction, the new EGRET limits on t he p ho- 
ton flux in the GeV region have severely limited [2a] the 
contribution of A-particle decay to the high energy end 
of the cosmic ray spectrum reported by the AGASA Col- 
laboration 29]. 

In our analysis we adopt the recent estimates of neu- 
tralino fluxes from decay of super heavy relics derived in 
Ref. [12J. The normalization is determined by matching 
the X-particle baryon flux to the difference between the 
observed spectrum at E « 10 11 GeV and contributions 
from a homogeneous population of astrophysical sources. 
Among the homogeneous models discussed in Ref. [l2j . 
the only one to accommodate both the ultra-high energy 
cosmic ray intensity and the GeV photon flux is a dis- 
tribution of X particles with rax ~ 2 x 10 12 GeV and 
primary 10-body decay X — > bq 5q. As can be seen in 
Fig. ^ this model is constrained to fit only the AGASA 
data below 10 11 GeV, with the baryonic flux from X- 
decay contributing less than 1/3 of the total [3(j ■ It is 
noteworthy that the associated flux of neutrinos from 
decay produced in this scenario, is also consistent with 
existing data. This is shown in Fig. [21 

As mentioned in the Introduction, in order to com- 
ply with existing fits of mSUGRA to the WMAP dark 
matter density, the contribution to this density from su- 
perheavy relics must be much less than that of the relic 
neutralinos. We recall that the cosmic ray flux result- 
ing from A-particle decays depends on the dimensionlcss 
parameter rx = £xto/i~ x , where £x = ^x/^cdm, and 
to is the age of the universe. Scenarios which include 
in the rate normalization photon flux from A-particles 
clustered in the halo lead to a value rx ~ 5 x 10~ n @. 
Omission of the photon channel in the normalization in- 
creases rx by a factor of about 10, and the extension 
from halo dominance to the homogeneous population 
used in this work further increases rx by a factor of 
15. Since models of X-production and decay typically 
lead to exponential dependence of both £x and t x on 
a reheating temperature Tr and a quantum mechanical 
tunneling action, respectively, there is no impediment 
on accommodating this change in rx while maintain- 
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FIG. 1: The solid line is a 2-component prediction of the ultra 
high energy proton flux, consisting of emission from ultra-high 
energy "stars" plus decay of super heavy relics, both sources 
distributed uniformly throughout the universe |30l) - The 
dashed line indicates the contribution from the X-decay, with 
initial state of 5 quarks + 5 squarks and mx = 2 x 10 12 GeV. 
The upper end of the spectrum as seen by different experi- 
ments (AGASA [llLHiRes [H, Fly's Eye [H SUGAR |H, 
and Haverah Park [3411 ') is also shown for comparison. 



ingf^x Ocdm ■ (For example, for the model provided 
in @, £x cx e - 2mx / TR ~ 1CT 4 - 1CT 8 .) 

Establishing a neutralino signal in upward going show- 
ers will require a suppression of neutrino events which 
create a background. This occurs naturally because the 
difference in neutrino and neutralino cross sections leads 
to differing absorption lengths in the Earth's crust. The 
neutrino flux is greatly attenuated by selecting events en- 
tering the Earth at angles > 5° below the horizon. For 
the EUSO effective aperture (duty cycle of 10%), the neu- 
trino background from a homogeneous population of X's 
decaying into bq bq (during the 3-yr mission lifetime) is 
about 0.3 events [33. The use of Poisson statistics im- 
plies the observation of > 3.09 events establishes a signal 
significance at 95% C.L. |3Sj . 

The determination of a neutralino event rate will de- 
pend strongly on the neutralino-nucleon cross section, 
and hence, through Eq. I|13|) . on the squark mass. Let us 
define P as the probability that the neutralino does not 
interact while traversing the Earth's crust. A large cross 
section will reduce P and lead to a lowering of emerging 
flux. In Fig.[3|we illustrate the behavior of P with energy, 
for two different values of m,q. On the other hand, it will 
enhance the event rate in the atmosphere for the neu- 
tralinos which do emerge without interacting. (Conser- 
vatively, we omit regeneration effects in the Earth when 
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FIG. 2: The solid (dashed) line indicates the neutralino [T^ ] 
(neutrino plus anti-neutrino |3C|) energy spectrum corre- 
sponding to the proton flux from X-particle decay given in 
Fig. Q The horizontally-hatched region at the top of the fig- 
ure has been already excluded from null results at AGASA + 
Fly's Eye (95% CL) [H and GLUE + FORTE (90% CL) [H . 
The thick dashed line indicates the 95% CL sensitivity of the 
EUSO mission. 



estimating a signal.) Apart from the direct effect on the 
size of the cross section, the squark mass will also influ- 
ence the event rate through the position of the resonant 
peak: a larger squark mass will probe higher value of 
E^ , and hence lead to a decrease in event rate because of 
the sharply falling flux. After considering these effects, 
we have found that a value of rriq ~ 1 TeV leads to an 
optimal signal. The degree of tolerance is fairly narrow: 
values lower than 800 GeV or larger than 1200 GeV will 
vitiate a signal at 95% CL. We note once more that this 
coincides with the region of Bino-dominance for the pa- 
rameter space compatible with WMAP data 0. 

For a given flux of neutralinos d$ / dE^ and observation 
time t, the event rate at EUSO is found to be 



N = 



dE^ N A P —a iN Ae DC t, (14) 



where Na is Avogadro's number and e DC « 10% is the 
duty cycle. The fraction of unscathed neutralinos, inte- 
grated over zenith angle < 85°, emerging upward in the 
direction of the orbiting telescope is found to be P m 
107 (E^/GeV)- - 256 . Now, by setting Ef a = 10 9 GeV 
and Ef 3 * = 10 11 GeV, one finds that for the full mis- 
sion lifetime EUSO will typically collect between 4 and 
5 events, and hence provide sensitivity to the neutralino 
flux given in Fig. [3 at more than 95% CL. 



V. CONCLUSIONS 
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FIG. 3: Fraction of neutralinos which traverse the Earth, for 
all zenith angles less than 85°, as a function of energy. 



Using a technique that exploits the different absorp- 
tion lengths of neutrinos and neutralinos in the Earth's 
crust, we have estimated the sensitivity of EUSO to 
isolate upward coming showers of ultra-high energy \- 
The neutralino-nucleon interaction has been approxi- 
mated by resonant squark production, with the neu- 
tralino being largely Bino in composition, and rriq ~ 1 
TeV. We have shown that, during the complete mission 
lifetime, the telescope will be sensitive to £"? dQ/dE^ > 

,-2 „„-i atthe 95% 



yr 



si- 



1.1 x 10~ 6 (Ex/GeV) - 7 GeV enr 
CL, for 10 9 GeV < < 10 11 GeV, and for the re- 



gion m q = 1.0 ± 0.2 TeV. A hard spectrum oc E- 1 - 3 
is typical of super heavy relic iV-body decays that are 
purely hadronic. This is a conservative estimate, since 
regeneration effects have been only considered in com- 
puting the neutrino background. We have explicitly an- 
alyzed the case in which the flux of ultra-high energy 
neutralinos is produced via decay of metastable heavy 
(mx = 2 x 10 12 GeV) particles with uniform distribution 
throughout the universe, and primary decay mode into 
5 quarks + 5 squarks. The normalization of tix/tx has 
been fixed to contribute about 1/3 of the events reported 
by the AGASA Collaboration below 10 11 GeV [ij. For 
this particular case, EUSO will collect between 4 and 5 
neutralino events (with 0.3 of background) in w 3 yr of 
running. 
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Existing limits on the diffuse photon flux in the GeV re- 
gion strongly limit the sensitivity of EUSO for primary 2- 
body decays of hadronic nature. This is because, for the 
normalization of the baryonic contribution to the ultra- 
high energy cosmic ray flux assumed in Ref. |12(, which 
is marginally consistent with new EGRET bounds, the 
accompanying neutralino flux produced for qq and qq is 
about an order of magnitude below the flux generated 
in the primary 10-body decay discussed above. On the 
other hand, the telescope can still be sensitive to the lep- 
tonic mode X — > 11. In this case, by reducing nx /tx by a 
factor of ~ 2, one lessens the problem with EGRET data 
and still leaves a window open for neutralino detection 
at EUSO. 

We note that for a Bino-like neutralino, the primary 
decay mode (90% branching fraction) of the squark is 
q — > qg, with a subsequent decay g — > qq \. Thus, the 
neutralino energy of the decay is about 1/6 of the primary 
energy. In the remaining 10% of the decays, q — ► qx- In 
either case, the shower energies are far above the ~ 1 
PeV threshold for the detector. If more detailed consid- 
erations are warranted in the future, regeneration effects 
during passage through Earth can be assessed, taking 
into account the energy losses of the decay modes. These 
effects will lead to some enhancement of P, and conse- 



quently of the event rate. 

We turn now to a brief discussion on the potential of 
the planned NASA mission Orbiting Wide-angle Light- 
collectors (OWL) [3^| . This mission will involve photo 
detectors mounted on 2 satellites in low equatorial orbit 
(600 - 1200 km). The eyes of the OWL will stereoscopi- 
cally image a geometric area of ~ 9 x 10 5 km 2 , yielding 
A ~ 3 x 10 6 km 2 sr. With its superior effective aper- 
ture (e DC « 10%), a 10 yr mission lifetime will allow one 
to discern on contributions of metastable relics to the 
upper end of the cosmic ray spectrum at the level of 1 
part in 10 2 . Consequently, the data from this mission will 
allow one to probe more deeply the parameter space of 
mSUGRA and its extensions. 
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